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Using first-principles calculations within density functional theory, we study the energetics and kinetics of 
carbon nucleation in the early stages of epitaxial graphene growth on three representative stepped metal surfaces: 
Ir(lll), Ru(OOOl), and Cu(lll). We find that on the flat surfaces of Ir(lll) and Ru(OOOl), two carbon atoms 
repel each other, while they prefer to form a dimer on Cu(ll 1). Moreover, the step edges on Ir and Ru surfaces 
cannot serve as effective trapping centers for single carbon adatoms, but can readily facilitate the formation of 
carbon dimers. These contrasting behaviors are attributed to the delicate competition between C-C bonding and 
C-metal bonding, and a simple generic principle is proposed to predict the nucleation sites of C adatoms on 
many other metal substrates with the C-metal bond strengths as the minimal inputs. 

PACS numbers: 68.35.Fx, 68.43.Bc, 68.43.Hn 



Since its first isolation, graphene has attracted rapidly grow- 
ing research interest because of its various intriguing prop- 
erties and potential applications in future electronics III |2|] . 
However, a route towards scalable mass production of qual- 
ity graphene for industrial use is still lacking. Among many 
newly developed techniques, epitaxial growth of graphene 
on metal surfaces offers a promising avenue |^ M HI B 
B S B E [ii E E H H E [l^ Large size and 
good quality graphene samples have been prepared on vari- 
ous metal surfaces The success in trans- 
ferring the epitaxial graphene grown on Ni and Cu surfaces 
to insulating substrates makes this method even more attrac- 
tive Hi. Additionally, various aspects about the growth 
mechanisms of graphene have been revealed in recent stud- 
ies of representative carbon/metal systems. For example, the 
growth of graphene on Ir(l 1 1) and Ru(OOOl) substrates is fed 
by the supersaturated two-dimensional (2D) gas of carbon 
adatoms, and a multi-carbon cluster attachment mechanism 
has been proposed with minimal effect of hydro- 
gen lITll [l3 1. On a Cu substrate, graphene is found to grow 
through a surface adsorption process, while on Ni it is by car- 
bon segregation or precipitation 

Despite these preliminary achievements, very little has been 
revealed about the growth kinetics, especially in the initial nu- 
cleation stages of carbon adatoms. Experimentally it has been 
found that carbon nucleation starts from the lower edges of 
steps on Ir(lll) [15] and Ru(OOOl) [10] surfaces, but it is un- 
clear why and to what extent nucleation at the step edges is 
preferred over terraces. Determination of nucleation sites is 
crucial in improving both the quality and quantity of epitax- 
ial graphene. In the growth of graphene on Ru(OOOl), mul- 
tiple nucleation on terraces can easily degrade the quality of 
graphene because defects will form at the interfaces of sep- 
arately nucleated graphene islands In graphene growth 
on Ir(l 1 1), the nucleation sites must not be too sparse, because 
otherwise rotated graphene domains are more likely to grow 
at the boundaries of the major phase of the islands that are 
aligned with the substrate |ll6l UJi . Quantity- wise, in order 



to eventually achieve mass production for industrial applica- 
tions, it is more desirable for nucleation of graphene islands 
to take place over the entire substrate rather than only at the 
edges of preexisting steps. In light of these aspects, a gen- 
eral guiding principle of determining the nucleation sites on 
different substrates will be highly beneficial. 

In this Letter, we present a comparative study of the ener- 
getics and kinetics in the initial stages of epitaxial graphene 
growth on three representative stepped metal surfaces, using 
first-principles calculations within density functional theory 
(DFT). We find that, whereas the interaction between two 
adatoms is attractive on flat Cu(l 1 1), leading to easy ad-dimer 
formation, it becomes repulsive on flat Ir(l 1 1) and Ru(OOOl), 
making ad-dimer formation improbable. On the other hand, 
even though the steps on Ir(lll) and Ru(OOOl) cannot serve 
as effective trapping centers for single carbon adatoms, such 
steps can readily facilitate the formation of carbon dimers at 
their lower edges. We rationalize these contrasting kinetic 
behaviors of carbon adatom diffusion and nucleation based 
on the delicate competition between the C-C bonding and C- 
metal bonding, and generalize this picture to predict the initial 
growth stages of graphene on different metal substrates. 

In our studies, we use the Vienna ab initio simulation pack- 
age (VASP) [18] with PAW potentials [19] and the general- 
ized gradient approximation (PBE-GGA) i20l] for exchange- 
correlation potential. All the metal surfaces are modeled by 
a 6-layer slab, with atoms in the lower 3 layers fixed in their 
respective bulk positions. A (2 x 4) squared surface unit cell 
is used to describe the Ir(lll) , Ru(OOOl) and Cu(lll) sur- 
faces. We use (322) and (332) surfaces to model the stepped 
Ir(l 1 1) and Cu(l 1 1) surfaces, which contains {100} (A-type) 
and {111} (B-type) microfacets, respectively. The stepped 
Ru(OOOl) surface is modeled by a vicinal surface with its nor- 
mal along the (0 11 10) direction, which contains alternating 
A- and B-type steps. All the terrace widths are ~1 1-12 A. The 
k-point mesh used in the calculations is 1 x 3 x 1 for stepped 
Ru(OOOl), and 3 X 3 X 1 for all the other cases [21]. We use the 
climbing image nudged elastic band (CINEB) method L22ll to 
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FIG. 1: (Color online) Top views of adsorption sites and binding 
energies of a C adatom around (a) Ir A-step, (b) Ir B-step , (c) Ru 
A-step, and (d) Ru B-step. The soHd curves represent C diffusion 
profiles. The vertical dashed line represents the position of step edge, 
and the horizontal dot-dashed line indicates the C binding energy on 
flat surfaces. Definition of labels: U - upper terrace, L - lower terrace, 
R - ridge site where C only binds to atoms in the ridge of a step edge; 
b - hep site, c - fee site. 



determine the energy barriers of the various kinetic processes. 

We first consider the adsorption and diffusion of isolated 
carbon atoms on flat metal surfaces. The most stable adsorp- 
tion sites and the corresponding binding energies, defined by 



C/subst 



Ec - Esuhsu on Ir(lll), Ru(OOOl) and 



hep (-7.66 eV) (in agreement 
l|]), and subsurface intersti- 



0, hep 



Cu(lll) are hep (-7.44 eV 
with previous calculations 
tial (-5.66 eV), respectively. The stronger binding on the other 
two substrates and the weaker binding on Cu(l 1 1) are consis- 
tent with the J-band model [23]: on Cu(lll), the C adatom 
mainly interacts with the free-electron like surface states of 
Cu, whose J-shell is completely filled; whereas on Ir(lll) or 
Ru(OOOl), the stronger binding originates from the hybridiza- 
tion between the sp orbitals of carbon and the half- filled ri- 
band of the substrate. On Ir(lll) and Ru(OOOl), the energy 
in the metastable fee sites are 0.25 eV and 0.74 eV higher, 
respectively. On Cu(lll), the metastable sites on the surface 
(fee, hep, bridge) are less stable than the subsurface interstitial 
sites by ~0.6 eV. The surface diffusion barriers (Sa) between 
a stable and the nearest metastable states are 0.75 eV, 0.87 eV 
and 0.66 eV on Ir(l 1 1), Ru(OOOl), and Cu(l 1 1), respectively. 

We next investigate the adsorption and diffusion of single 
C adatoms at step edges of Ir(l 1 1) and Ru(OOOl). The results 
for the case of Cu(lll) will be reported in a separate work 
1^4!], as will be explained later. As shown in Fig.[Tl the calcu- 
lated binding energies at step edges are not much larger than 
those on flat surfaces. The same is true for the kinetic barri- 
ers. Considering the high growth temperatures in experiments 
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FIG. 2: (Color online) Binding energies of two C adatoms on flat 
metal surfaces as a function of their separation distance. Data points 
around the vertical dashed line correspond to the formation of C 
dimers. Inset shows the top view of a C dimer on a close-packed 
metal surface. Kinetic barriers are not shown. 



(~1000 K), we arrive at the conclusion that the substrate steps 
do not serve as effective traps for single C adatoms. The ab- 
sence of large step-crossing barriers and deep wells at step 
edges, in contrast to the traditional Ehrlich-Schwoebel (ES) 
picture |[25l|2^[27|], is attributed to the passivation of the least 
coordinated step edge atoms by the non-metal adatoms, lead- 
ing to a large amount of energy gain ll28l] . Therefore, in the 
transition state where the adatom bonds to step edge atoms 
only, it more effectively passivates the edge atoms, resulting 
in a lower barrier. Similarly, because at lower step edges the 
substrate atoms are relatively over-bonded, C adatoms cannot 
gain much energy by having more neighbors 1,28,1 . Thus no 
extra deep potential well is present either. 

Where should nucleation of carbon adatoms occur if they 
are not effectively trapped anywhere on the substrates of Ir 
and Ru? The above results indicate that knowing the behav- 
ior of non-interacting single carbon adatoms is insufficient 
to answer this question. Specifically, since carbon adatoms 
are known to form strong covalent bonds with one another 
when they nucleate to form graphene, it is necessary to take 
the carbon-carbon interaction into account. We therefore next 
study the formation of C dimers as the first step of nucleation 
on the metal substrates. Fig. [2] shows the trend of binding 
energies of two carbon adatoms on the flat metal surfaces. 



defined by A^2C = E'^ 



2C/subst 



2Ec - £^subst, as a func- 



tion of the separation distance. One can immediately notice 
that on Ir(lll) and Ru(OOOl) the formation of C dimers is 
energetically unfavorable, but on Cu(lll), dimers are much 
more stable than separate C adatoms by over 2 eV. More- 
over, the energy barrier of forming a dimer for two neigh- 
boring C adatoms is only 0.32 eV on Cu(l 1 1), which is much 
smaller than those on Ir(lll) (1.37 eV) and Ru(OOOl) (1.49 
eV). These findings suggest that on Ir(lll) and Ru(OOOl), 
C adatoms are mutually repulsive and cannot form dimers, 
whereas on Cu(lll) they strongly attract each other, leading 
to the formation of dimers and larger islands. 
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FIG. 3: (Color online) (a) and (b) Top- and side-view of the most 
stable configuration of a C dimer at the lower edge of a (a) A-step and 
(b) B-step. (c-h) Binding energies of two carbon adatoms with one 
C atom fixed at the lower step edge and another moving on the upper 
terrace, the lower terrace and along the lower step edge, respectively. 
Horizontal axis is their separation distance. Vertical dashed line in 
each panel shows where a C dimer is formed, and horizontal dot- 
dashed line shows the binding energy of two separate C adatoms on 
flat surfaces. 



TABLE I: Binding energy difference between a C dimer at lower step 
edges (^dimer/step) ^nd (1) two separate C adatoms on flat surfaces 
(^2C/flat) and (2) a C dimer on flat surfaces (^dimer/flat). (The unit 
is in eV). 



Energy difference 


Ir 




Ru 




Cu 




Step 


A 


B 


A 


B 


A 


B 


-^dimer/step -^dimer /flat 


-1.31 - 


1.19 


-1.09 - 


0.68 


-0.64 - 


0.28 


-^dimer/step ~ -^2C/flat 


-0.44 - 


0.33 


-0.46 - 


0.15 


-2.54 - 


2.18 



Now that the nucleation sites on Cu(lll) have been iden- 
tified, we next show that on Ir(lll) and Ru(OOOl) nucleation 
can be readily facilitated by the step edges. Table H compares 
the binding energies of a C dimer at lower step edges with the 
cases of a dimer on flat surfaces and two separate C adatoms 
on flat surfaces, showing that dimers at step edges are not only 
much more stable than on flat surfaces, but also more stable 
than two separate C adatoms. Therefore, even though on flat 
Ir and Ru surfaces C dimerization is not preferred, C adatoms 
attract each other at lower edges of the surface steps. To better 
illustrate this, we plot the binding energies of two C adatoms 



on stepped metal surfaces with their separation in Fig. O In 
all cases, there is a deep potential well upon the formation of 
a C dimer at lower step edges. 

Summarizing the above results, we have shown that on 
Ir(lll) and Ru(OOOl), nucleation of C adatoms first occurs 
at substrate step edges, in agreement with existing experi- 
ments ll5|]; whereas on Cu(lll), our results predict that 
C adatoms should nucleate everywhere on the surface. We 
note that even though there is also a deeper potential well for 
the C dimer formation at the Cu(l 1 1) steps, such steps are not 
so crucial in the nucleation of C adatoms on Cu, because C 
adatoms are already strongly attractive to each other on the 
terraces and readily form dimers before they can reach a step 
edge. For the same reason we ignored the discussion about 
Cu(l 1 1) steps in the study of single C adatoms above. 

For C adatoms on Ir(lll) and Ru(OOOl), this exceptional 
tendency towards dimerization at substrate step edges is re- 
lated to the special local bonding geometry of a C dimer at 
those sites. In Figs. [3] (a) and (b) the bonding geometries of 
a C dimer at A-type and B-type step edges are shown, re- 
spectively. By comparing those with a C dimer on flat sur- 
faces shown in the inset of Fig. [2l one can observe that the 
bonds in the latter case are severely twisted. Since the cova- 
lent bonds are highly directional and it is energetically costly 
to change the relative bond angles, the relaxation of the cova- 
lent bonds by the step geometry leads to the extra stability of 
the C dimers. 

The contrasting behavior of the interacting C adatoms on 
flat close-packed Ir(lll), Ru(OOOl), and Cu(lll) surfaces can 
be attributed to the competition between the C-C and C-metal 
interactions. The C-C bond lengths of carbon dimers on flat 
Ir(lll), Ru(OOOl) and Cu(lll) surfaces are 1.397 A, 1.376 
A, and 1.299 A, respectively, which are very close to the 
length of a C-C double bond (1.34 A). A double bond re- 
quires two bonding electrons from each C adatom, but one 
carbon adatom has only four valence electrons and three near- 
est metal neighbors on the surface. So intuitively, the forma- 
tion of a C dimer will weaken the C-metal bonding because 
of less bonding electrons. Therefore, if the C-metal bonds are 
very strong, which is the case of Ir and Ru, the dimer forma- 
tion is not energetically favorable. Conversely, in the case of 
Cu where C-metal bonding is weak, formation of a dimer is 
preferred for two C adatoms. 

Next we show that the above picture is not limited to the 
three representative cases, but can be generalized into a simple 
guiding principle. To this end, we compare the binding ener- 
gies of C adatoms and C dimers on the close-packed surfaces 
of various transition metals, as shown in Fig.Hl It is apparent 
that the weaker the C-metal interaction is, the more preferred 
the C dimers are. In all the cases of noble metals, which have 
closed J-shells and strong free-electron like surface states, C 
dimerization is preferred. The relative strength of C-metal 
and C-C interactions largely determines whether the net in- 
teraction between C adatoms is attractive or repulsive. The 
dimer-preferred and dimer-not-preferred systems are essen- 
tially separated by the vertical dashed line corresponding to 
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FIG. 4: (Color online). Binding energy difference between a C 
dimer (A^^dimer) and two separate C adatoms (2 AEc) with respect 
to AEc on close-packed transition metal surfaces. Vertical dashed 
line corresponds to the binding energy of a C-C double bond (-6.33 
eV)L2a]. 



the energy of a C-C double bond (-6.33 eV) |29]. The de- 
viation from this trend may be, for example, because of the 
variation in bonding nature or geometrical effects. Based on 
the results presented earlier and the prototypical nature of the 
systems we have studied therein, we can further conclude that 
for those systems in which C dimer s are not preferred on ter- 
races, C nucleation should first occur at substrate step edges. 
Thus, our study makes it possible to predict where the initial 
nucleation should happen armed solely with the knowledge of 
the binding energy of C adatoms to the metal substrate. 

This generic principle can lead to many strong predictions. 
For example, in the strong C-metal binding regime, a flat sub- 
strate with scarce steps may not result in growth of quality 
graphene because of the simultaneous nucleation at multiple 
sites on the terraces [10], a somewhat counterintuitive conclu- 
sion. In the weak C-metal binding regime, epitaxy on single- 
crystal flat Cu(lll) is more likely to yield graphene with the 
desired high quality and potential mass production, because C 
adatoms prefer to nucleate everywhere, and the mismatch of 
graphene with Cu substrate is very small. 

We finally emphasize on another salient prediction of the 
present study. The most stable configuration of a C dimer 
at A- or B-type step edges are shown in FigsIS (a) and (b), 
with the dimers bridging the step ridge and the lower terrace 
perpendicular to the step. Experimentally it is already possi- 
ble to prepare near-perfect vicinal surfaces containing regu- 
lar straight step arrays by precise control of the miscut [30]. 
Therefore on Ir(lll), Ru(OOOl), and other metal surfaces with 
strong C-metal bonding, growth of straight graphene nanorib- 
bons with zigzag edges are likely to be obtained from the 
regularly spaced C dimer arrays formed along the A- or B- 
type steps in the nucleation stage. This suggestion may of- 
fer a potentially more attractive route for mass production of 
highly ordered graphene nanoribbons with zigzag edges than 



the known approaches ||31l |32 

In summary, we have performed a comparative study of the 
energetics and kinetics of carbon adatoms on stepped Ir(lll), 
Ru(OOOl) and Cu(lll) surfaces, with intriguing predictions. 
We have found that the substrate step edges cannot effectively 
capture single C adatoms, and two carbon adatoms dislike 
forming dimers on flat Ir(lll) and Ru(OOOl) surfaces either, 
though they strongly attract each other on Cu(lll). How- 
ever, the lower edges of steps on Ir(lll) and Ru(OOOl) can 
readily mediate the formation of C dimers. These findings 
have been rationalized by considering the competition be- 
tween the carbon-carbon and carbon-metal bonds. We have 
further generalized this picture to a simple guiding princi- 
ple for predicting the initial nucleation sites of C adatoms on 
many other metal substrates, and the predicted behaviors in 
the early stages of epitaxial graphene growth are expected to 
be instrumental in achieving mass production of high quality 
epitaxial graphene. 
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